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PHYSICAL NATURE OF HYDROGEN BOND 


The physical nature and the correct definition of hydrogen bond (H-bond) are considered. The 
influence of H-bonds on the thermodynamic, kinetic, and spectroscopic properties of water is 
analyzed. The conventional model of H-bonds as sharply directed and saturated bridges between 
water molecules is incompatible with the behavior of the specific volume, evaporation heat, and 
self-diffusion and kinematic shear viscosity coefficients of water. On the other hand, it is shown 
that the variation of the dipole moment of a water molecule and the frequency shift of valence 
vibrations of a hydroxyl group can be totally explained in the framework of the electrostatic 
model of H-bond. At the same time, the temperature dependences of the heat capacity of water 
in the liquid and vapor states clearly testify to the existence of weak H-bonds. The analysis 
of a water dimer shows that the contribution of weak H-bonds to its ground state energy 
is approximately 4~5 times lower in comparison with the energy of electrostatic interaction 
between water molecules. A conclusion is made that H-bonds have the same nature in all other 
cases where they occur. 

Keywords: hydrogen valence vibrations of a water molecule, frequency shift, hydrogen bond, 
electrostatic origin. 


The hydrogen bond is a phlogiston of the 20-th century. 

L. BULAVIN 

1. Introduction 

During the last 100 years, the concept of hydrogen 
bond played a large role in the description of proper¬ 
ties of water, alcohols, aqueous alcohol solutions, and 
other systems, in which hydrogen bonds were sup¬ 
posed to exist. Unfortunately, a consistent theory for 
the latter has not been created within this time in¬ 
terval. 

Hydrogen bonds obviously neither belong to stan¬ 
dard chemical bonds with energies of about lOOfeeTtr, 
where /cb is the Boltzmann constant and Ttr is the 
temperature of water’s triple point, nor to ionic ones 
with energies of the same order of magnitude. At the 
same time, the energies, which are ascribed to hydro¬ 
gen bonds, exceed those of van der Waals or disper¬ 
sion forces by 1-2 orders of magnitude. 

M.D. Sokolov was the first who paid attention to 
this circumstance almost 60 years ago [l|. He indi¬ 
cated that the energy of hydrogen bonds approaches 
the energy of dipole-dipole interaction. In such a way. 


M.D.Sokolov showed that the main contribution to 
the hydrogen bond energy is made by the electrostatic 
interaction, which does not associated with the de¬ 
formation of electron shells in a water molecule. This 
idea was also supported by an outstanding physical 
chemist Prof. G.G. Malenkov during oral discussions 
on the hydr og en bond problem. Later, it was shown 
in works [2 t 2| that the deformation of electron shells 
in water molecules forming a hydrogen bond is a small 
quantity and generates an irreducible contribution to 
the energy of interaction between molecules. It is this 
contribution that should be called the hydrogen bond. 

In the work by Clement! Q , the optimum configu¬ 
rations of hydrogen and oxygen atoms in the isolated 
water molecule and in the water molecule located near 
the Li’*' ion were calculated. The difference between 
the oxygen-hydrogen distances in those two configu¬ 
rations amounts to about 0.005 A. This fact clearly 
testifies in favor of water molecule models with fixed 
invariant positions of model charges, because the dis¬ 
tances between the centers of charges in the water 
molecule did not change even under the action of an 
electric field with a high strength that occurred near 
the Li+ ion. 
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The facts indicated above substantiate the appli¬ 
cation of a number of model electrostatic potentials 
proposed or described in works Those po¬ 

tentials allow one to explain the thermodynamic and 
the majority of kinetic properties of water, as well as 
many other systems. At the same time, the temper¬ 
ature dependence of the heat capacity of water 171 
undoubtedly testifies to the necessity to consider the 
contributions made by hydrogen bonds (the energy of 
these bonds is about (1.5-^3)A:B7tr). A new view on 
the character of the intermolecular interaction also 
follows from the systematization of liquids proposed 
by Leonid Bulavin (see also works diml). 

In this work, we tried to systematize all main facts 
that testify to the necessity of a radical revision of 
the concept of hydrogen bonds. Our conclusions are 
based on the results obtained for the temperature de¬ 
pendences of the self-diffusion and kinematic shear 

viscosity coefficients, specific volume j _ evaporation 

heat [ 23 , and heat capacity of water 1171.12.31.12411. In 
addition, we analyze the frequency shift for valence 
hydrogens in the water molecule |25l | and the variation 
of its dipole moment under the action of environmen¬ 
tal molecules Q. 


2. Manifestation of Hydrogen 
Bonds in Transport Processes 

In this section, the manifestations of hydrogen bonds 
in two most important kinetic transport phenomena 
are discussed. These are self-diffusion and shear vis¬ 
cosity. Note that, in order to exclude the density ef¬ 
fects, the behavior of kinematic shear viscosity will 
be analyzed. The main results are obtained with the 
use of experimental data measured at the coexistence 
curve. 

From the conventional viewpoint, hydrogen bonds, 
which couple neighbor molecules with an energy of 
about lOfeeTtr, should interfere with the thermally 
driven drift of molecules and the relative shift 
of adjacent liquid layers. On the other hand, the 
rotational motion of molecules has to be taken into 
consideration, which follows from the basic principles 
of statistical mechanics. It should be noted at once 
that those two statements do not agree with each 
other. If hydrogen bonds are imagined as rods or 
string segments, the rotational motion of molecules 
becomes impossible. In this connection, an assump¬ 
tion should be made that the transverse elasticity of 
hydrogen bonds either equals zero or its value does 





Fig. 1. Dependence of the ratio between the dif¬ 

fusion coefficients of water and argon at the corresponding tem¬ 
peratures 


not exceed feeTtr. The concept of hydrogen bond as 
a physical object with so controversial properties is 
unsatisfactory. 

In order to confirm a conclusion about an unsat¬ 
isfactory character of such notions concerning the 
hydrogen bond, let us consider the temperature de¬ 
pendences of the relative values of the self-diffusi¬ 
on and viscosity coefficients for water (w) and ar¬ 
gon (Ar). The temperature behaviors of the self-dif¬ 
fusion, and kinematic shear viscosity, 

£'(“')coefficient ratios are depicted in Figs. I 
and 2, respectively (here, A = A/Atr). Attention 
should be paid that the coefficients are compared at 
the same relative temperatures t = T/T^l\ where 
is the corresponding triple point temperature of 
the liquid {i = w or Ar). The states of argon and 
water in the interval from the triple point to the crit¬ 
ical one and supercooled water are considered. The 
values of self-diffusion and viscosity coefficients for 
water were taken from the NIST reference database 
(26| . and the self-diffusion coefficients for argon were 
calculated, by using the molecular dynamics methods 
because of the lack of detailed experimental 


27, 


results. 

As one can see from the dependences exhibited in 
Figs. I and 2, hydrogen bonds practically do not 
manifest themselves in the temperature behavior of 
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I- 



ficients of kinematic shear viscosity in water and argon at the 
corresponding temperatures 

the diffusion and viscosity coefficients. First, the 
self-diffusion coefficient of water exceeds that of ar¬ 
gon, which completely agrees with a higher mobility 
of lighter water molecules (by order of magnitude, 
the self-diffusion coefficient Dg k-QT/v). Second, 
the viscosity coefficient for water is practically by 
an order of magnitude smaller than that of argon, 
which also contradicts the initial statement that hy¬ 
drogen bonds may strongly affect the viscosity of wa¬ 
ter. Moreover, in the whole temperature intervals of 
the liquid state existence, the self-diffusion and vis¬ 
cosity coefficients of both water and argon satisfy the 


Table 1. Self-diffusion coefficient 
of water molecules in aqueous so lutions 
of single-charged electrolytes fe 





LiBr(30) 

Lil(24.8) 

X 10®, cm^/s 



2.3 

2.18 



NaCl(15.9) 

NaBr(16.5) 

Nal(16.1) 

X 10®, cm^/s 


2.14 

2.26 

2.38 


KF(15) 

KC1(16.1) 

KBr(16.1) 

KI(16.4) 

X 10®, cm2/s 

1.99 

2.44 

2.68 

2.8 


relation 

- l/r„, (1) 

where is the radius of the molecule, i.e., this com¬ 
bination is close to the corresponding constant for 
argon or water. The invariance of this combination 
for water testifies that hydrogen bonds have no sub¬ 
stantial relation to the problem of the water viscosity 
behavior. 

This conclusion evidently correlates with the state¬ 
ment made in work that the behavior of viscosity 
is governed by the averaged interaction potential be¬ 
tween molecules. The averaging is a consequence of 
the almost free rotation of water molecules, which 
would be impossible in the case of rod-like hydrogen 
bonds. 

Moreover, we attract attention to the absence of a 
network formed by hydrogen bonds. This conclusion 
immediately results from a careful analysis of the self¬ 
diffusion coefficients and the mobilities of ions and 
water molecules in diluted electrolyte solutions (^ . 

3. Mobility of Water 

Molecules in Electrolyte Solutions 

In this section, the physical nature of the mobility 
of water molecules in diluted aqueous solutions of 
electrolytes, when there are no more than 15 water 
molecules per ion, is discussed. Attention is focused 
on the fact that the behavior of the mobility coef¬ 
ficients of water molecules-or, in other words, their 
self-diffusion coefficient-is ultimately determined by 
the radii of ionic hard cores. Hence, the network of 
hydrogen bonds, the existence of which is postulated 
in the overwhelming majority of works, does not man¬ 
ifest itself. 

3.1. Self-diffusion coefficients 
of water molecules 

First of all, we would like to attract attention to how 
the behavior of the self-diffusion coefficients of water 
molecules depends on the dimensions of cations and 
anions. The corresponding values of self-diffusion co¬ 
efficients of water molecules in several diluted 

electrolyte solutions at the temperature T = 296 K 
are listed in Table 1. In each of three table rows, the 
cation remains the same, i.e.. Table I exhibits the 
dependence of on the anion size. The solution 
concentration in the table is presented by the number 
Zw of water molecules per ion. The Zw-yahie is shown 
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in parentheses near the chemical formula of elec¬ 
trolyte, e.g., NaCl(15,9). Hence, the minimum value 
Zuj = 15 corresponds to an electrolyte concentration 
of 3.3 moL%. At such a concentration, the mutual in¬ 
fluence of cations and anions can be neglected with a 
satisfactory accuracy. Table 1 should also include the 
self-diffusion coefficient of water molecules in the so¬ 
lution Csl(17,4), namel y, = 2.71 x 10® cm^/s. In 
accordance with works 3^l33| , the self-diffusion coef¬ 
ficient of molecules in water at the same temperature 
T = 296 K equals 

= 2.35 X 10“® cmVs. (2) 


One can see that the self-diffusion coefficients of 
water molecules grow as the sizes of cations (the data 
in columns) and anions (the data in rows) increase, 
which could testify to a destruction of the hydrogen 
bond network. At the same time, in this case, water 
molecules would have freely surrounded electrolyte 
ions to form hydration spheres and to give rise to a 
substantial reduction of the self-diffusion coefficients 
of water molecules, which contradicts the data quoted 
in Table 1. 


3.2. Features in the dependence 

of the self-diffusion coefficients of water 

molecules on the ion size 


To make the analysis of the results quoted in Ta¬ 
ble 1 more comprehensive, let us consider the corre¬ 
lations between the self-diffusion coefficients of wa¬ 
ter molecules in electrolyte solutions and the radii of 
dissolved ions. The considered radii of ions (i) were 
determined from crystallographic data, (ii) were se¬ 
lected in a way to favor the correct reconstruc¬ 
tion of the molecular dynamics in electrolyte solu¬ 
tions with the use of computer simulation methods, 
and (iii) were evaluated from the ionic polarizability 
values. 

The first row in Table 2 contains the values of crys¬ 
tallographic radius Vc of ions [^. The second row 
contains the radii of ions determined in computer 
experiments aimed at the description of the disper¬ 
sion (van der Waals) interaction between ions and 
water molecules • The radii of ions , which were 
determined from ionic polarizabilities a, by using the 
formula [s^ 

(3) 

are quoted in the third row. The radii Vc, r„, and 
will be referred to as hard ionic radii. 


From Tables 1 and 2, it follows that the de¬ 
pendences of the self-diffusion coefficients of water 
molecules on the radii of ionic hard cores demonstrate 
the following regularities: 

1 ) for diluted lithium and sodium electrolyte so¬ 

lutions, in which Vf. < ^1^, where is the aver¬ 
age distance between the oxygen atoms in neighbor 
water molecules, the inequality Ds’“^(el) < is 

obeyed. The inequality Tc < \lw is obviously vio¬ 
lated only for Cs’*'. In the solutions of potassium elec¬ 
trolytes, in which rc(K+) ~ \lw^ a transition from the 
previous inequality between the self-diffusion coeffi¬ 
cients of water molecules to the inequality (el) > 
> is observed; 

2 ) in electrolyte solutions with a fixed cation, but 
for the lithium one, the self-diffusion coefficients of 
water molecules grow together with the anion radius; 

3) in the lithium electrolytes, the character of the 
dependence of Ds’“^(el) on the anion radius is opposite 
to that described in the previous item. 

Since the concentrations of different single-charged 
cations and anions are close to one another, the dif¬ 
ferences in the behavior of Ds'^^(el) can be associ¬ 
ated with geometrical factors and their different ac¬ 
tion on the structure of local environment (hydration 
effects). The former possibility should be rejected, be¬ 
cause the geometrical obstacles should diminish with 
the reduction of the cation radius, which contradicts 
experimental data. At the same time, a local restruc¬ 
turing of water in close vicinities of cations and anions 
turns out more appreciable for larger radii of their 
hard cores. 

It should be noted that the effect of local restruc¬ 
turing in the solution is insignificant, because an in¬ 
crease or a reduction of the self-diffusion coefficient 
for water molecules in most cases does not exceed 
10 % and is proportional to the molar concentration of 


Table 2. Hard core radii (three upper 
rows) and Stokes radii of cations and anions 



Li+ 

Na+ 

K+ 

Cs+ 

F- 

Cl- 

Br- 

I- 

rc, A 

0.6 

0.95 

1.33 

1.69 

1.36 

1.81 

1.95 

2.16 

Ta, A 

0.76 

1.3 

1.67 

1.94 

1.56 

2.2 

2.27 

2.59 

Ta, A 

0.45 

1.12 

1.41 

2.02 

1.51 

2.33 

2.55 

2.93 

A 

2.38 

1.84 

1.25 

1.19 

1.66 

1.21 

1.18 

1.19 

A 

1.91 

1.88 

1.14 

1.15 

1.77 

1.33 

1.22 

1.35 
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electrolyte admixtures. All those facts are completely 
incompatible with the statement about the presence 
of a developed network of hydrogen bonds in water 
and diluted aqueous electrolyte solutions. 


3.3. Nonideality degree of electrolyte solutions 

The results obtained should be appended with the 
data concerning the nonideality degree of diluted 
electrolyte solutions, which were obtained in work 
[ 2 ^ . According to the cited work, the nonideality de¬ 
gree can be calculated using the relation 


(5 



X 





( 4 ) 


where pw'^ and are the densities of water and 
electrolyte, respectively, in the liquid or amorphous 
state, and p is the density of the aqueous electrolyte 
solution. The densities of some electrolyte solutions 
can be found in Table 3, and the degrees of their 
nonideality, 6, are given in Table 4. 

Note that the values of nonideality parameter ob¬ 
tained by formula ([4]) are averaged over the number 
of water molecules that surround the corresponding 
ion. One can see that relatively small values of the 
mobility of lithium cations and the self-diffusion co¬ 
efficient of water molecules in the lithium electrolytes 
correspond to negative nonideality degrees. This fact 
testifies that lithium cations do not promote the for¬ 
mation of hydration spheres around them with the 


Table 3. Densities of aqueous electrolyte 
solutions at the fixed concentration aigi = 4 wt.% 



F- 

cr 

Br- 

I- 

Li+ 


1.02 

1.031 

1.033 

Na+ 


1.026 


1.035 

K+ 

1.032 


1.029 

1.032 

Cs+ 

1.034 

1.035 




Table 4- Parameters of solution 

nonideality at the fixed concentration x^i = 4 wt.% 



cr 

Br- 

r 

Li+ 

-0.001 



Na+ 

0.003 



K+ 


0.003 

0.004 

Cs+ 

0.005 




density exceeding that of water. In other cases, if one 
may talk about the formation of hydration spheres, 
the smallness of the parameter S testifies that their in¬ 
fluence on the density in diluted electrolyte solutions 
is weak. 

From the facts mentioned above, it follows that 

(i) the key role in governing the transport prop¬ 
erties of aqueous electrolyte solutions-first of all, the 
behavior of the mobility coefficients of ions and water 
molecules-is played by hard cores of those objects; 

(ii) the conventional scenario that ions in aqueous 
solutions move through “voids” existing in the hydro¬ 
gen bond network is incorrect; 

(iii) hydrogen bonds between water molecules have 
to be considered as a convenient model, which to a 
certain extent reflects the existence of correlations be¬ 
tween the dipole moments of molecules, as well as 
multipole moments of higher orders; and 

(iv) the role of hydration effects is insignificant and 
can be taken into account in the framework of the 
thermodynamic perturbation theory. 

4. Experimental Evidence 

for the Similarity of the Thermodynamic 

Properties of Water and Argon 

In this section, some facts testifying to the thermo¬ 
dynamic similarity between water and argon and, 
hence, calling into question the conventional view¬ 
point about the crucial role of hydrogen bonds in 
the formation of water properties are presented. With 
that end in view, let us consider the temperature de¬ 
pendences for the simplest quantities: the fractional 
volume (this is one of the mechanical characteristics 
of the system) and the evaporation heat (the most im¬ 
portant among the thermal parameters). In parallel, 
the behavior of ordinary water and its heavy coun¬ 
terpart, which considerably differ from each other by 
the character of molecular rotational motion, will be 
analyzed. 

The temperature dependences of the fractional vol¬ 
umes of water {i = w) and argon (i = Ar) 
at their coexistence curves will be compared in the 
spirit of the similarity principle for the corresponding 
states of the system [^. This means that the ratio 
j^(H 20 )(^) — y{ii 20 ) (j.yyiAr) between the normal¬ 
ized volumes where means 

the (t)-value at the corresponding critical point, 
should be considered as a function of the dimension- 
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less temperature t = T/Tc''\ where t = T/Tc'"'^ is the 
critical temperature of the i-th liquid. As is seen from 
Fig. 3, practically in the whole temperature interval of 
water existence in the liquid state, 0.42 < t < 0.9, the 
temperature dependences of the fractional volumes of 
water and argon are similar. Their ratio in 

the temperature interval 0.55 < t < 0.8 is approxi¬ 
mated with a satisfactory accuracy by the linear de¬ 
pendence 

= (5) 

with the coefficients 

H 2 O : = 0.827 = 0.021 = 0.83, 

( 6 ) 

D 2 O : = 0.851 = 0.021 = 0.83, 

i.e. it remains almost constant. 

The deviations of the ratio 
/r;(H 20 )(^) from unity is appreciable only for at t 
<C 0.55 and in a vicinity of the critical point. How¬ 
ever, they do not exceed 4%. 

A comparison between the evaporation heats 
for water and argon is even more intriguing (see 
Fig. 4). The deviation of from an approxi¬ 

mate value of about 6.2 does not exceed 1.1%. 

Analogously to formula ([5]), the ratio (t) in 

the temperature interval 0.55 < t < 0.8 is quasilinear. 




( 7 ) 


where 

H 2 O : Af'’ = 6.137 = 0.080 A, = 0.85 

D 2 O : rlf'> = 6.084 Af^ = 0.182 A„ = 0.85. 


( 8 ) 


Moreover, the linear functions 1 — and 1 — Xqt 
in Eqs. (jS]) and (O, respectively, are almost identical, 
which testifies to their common origin. In works |17l 
23|, it was shown that those functions are generated 
by weak hydrogen bonds. 

Small deviations of Rv\t) and R^q\t) from their 
constant values evidently testify to a weak effect of 
hydrogen bonds and a similarity of intermolecular 
potentials in water and argon. The latter conclusion 
has a completely natural explanation. The behavior 
of the fractional volume and the evaporation heat for 
water is governed by the averaged potential of inter¬ 
action between molecules. In its turn, the potential 



Fig. 3. Temperature dependences of the ratio = 

= calculated at the coexistence curves of 

water and argon in accordance with experimental data l26ll 



Fig. 4- Temperature dependences of the ratio = 

= q(H 20 ) !q(A-T) accordance with experimental data l2fill 


averaging is a direct consequence of the rotational mo¬ 
tion of water molecules. As a result, the anisotropic 
effects associated with weak hydrogen bonds become 
practically smoothed out. 

5. Vibration Frequencies 

of Hydrogen Atoms in the Water Molecule 

in Vapor, Water, and Ice 

Let us discuss the frequency shift for the longitudinal 
(valence) vibrations of the hydrogen ion in a wa¬ 
ter molecule, which lies close to the line connecting 
the centers of mass of oxygen atoms in two neighbor 
water molecules composing a dimer (see Fig. 5). The 
corresponding vibration frequencies are determined 
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and the reduced mass of the oxygen-hydrogen system 
in the water molecule approximately equals 



Fig. 5. Configuration of a linear dimer of water molecules. 
The equilibrium angle values are y = 3.76° and 9 = 41.1° 



Fig. 6. Dependences of the relative change AK^rV^rr^ of the 
force constant on the distance roi 02 the standard dimer: 
with regard for a mutual tuning of molecular orientations (1, 
triangles) and for the fixed orientation of molecules as in the 
configuration in Fig. 5 (^, curve) 

by the formula 


^11 



( 9 ) 


Mrb 


(Mo + Mh) Mh 
(M o + Mh) + Mh 


( 11 ) 


The approximate character of Eq. m is explained 
by the fact that hydrogen (see Fig. 5) is not lo¬ 
cated at the line connecting the centers of mass of 
oxygen atoms. The double derivative in Eq. (HUD is 
assumed to be taken at the point determined by 
the equation 




$ 


f^(rHi,rH2) + $Int(fOiH; 


3 ? 1*OiH4 5 



5fHi 


= 0 . 
( 12 ) 


It should be noted that Eq. m gives rise to only 
an insignificant displacement Achi of the equilibrium 
position of hydrogen . Assuming that 


Ihi =fHi + AfHi, 
where 


|AfHi| < fHi, 

the value of AJhi can be determined with the help of 
a more simple equation, 


K^jAfHi + Vfjj^<i)int(fo 1H3J f'OiH4: i’0i02) 


= 0 , 


( 13 ) 


where the force constant (the dimensionless quantities 
$ = ^/kBTtr and where a = 298 A is 

the oxygen diameter, will be used) is defined in a 
standard way: 



0=2 I 


( 10 ) 


where is the elasticity coefficient for the bond 
between the hydrogen and oxygen atoms in the 
monomer. 

One may get convinced (see work 0) that, by 
the order of magnitude in the framework of the elec¬ 
trostatic model, the ratio Afni/fni, where AJhi = 
= — thi ) satisfies the inequality 


Table 5. Force constants for water molecules 





a^s 


drf 


dridO 

GSD 

256.98 

182.15 

33.439 

Experim. 

256.98 

190.64 

33.439 


AfHi/fHi (fHi/fOiOa)"^ < 0.02. 

Therefore, the displacement of a hydrogen atom can 
be neglected practically for all distances between the 
oxygen atoms in the water dimer. 

The values of force constants for the water 
molecule, which correspond to the GSD potential, are 
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listed in Table 5. For the sake of comparison, Table 5 
also contains the corresponding force constants deter¬ 
mined experimentally. 

The force constant for symmetric valence vibra¬ 
tions in the equilibrium configuration of a dimer 
shown in Fig. 5 can be calculated by formula ([TUI) . 
The dependence of the relative change of a force con¬ 
stant on the distance between the oxygens in the 
dimer is depicted in Fig. 6. Point W with the coor¬ 
dinates (0.939, —0.051) corresponds to the distance 
roi 02 = 2.8 A between the oxygen atoms, which 
is typical of water near its triple point. The rela¬ 
tive change of the constant of symmetric valence vi¬ 
brations at this point equals AK^rV^rr^ = —0.05. 
Hence, the values of the force constant of valence vi¬ 
brations at the distances between the oxygen atoms 
corresponding to a dimer at equilibrium and to liquid 
water differ from each other by 5%. 

5.1. Results of calculations 

of the valence vibration frequencies 

for hydrogen atoms in the water molecule 

The frequency shift of hydrogen valence vibrations 
in the water molecule depends on the water phase 
state and can reach several hundreds of inverse cen¬ 
timeters (Table 6). In this work, we assume that the 
main contribution to the experimentally observed fre¬ 
quency shift is made by electrostatic forces associated 
with multipole moments of water molecules. 

The major result of our research consists in that the 
electrostatic forces really induce the frequency shifts, 
which agree with experimental data by both the shift 
direction and the order of magnitude. The frequency 
shift of valence vibrations equals 

. 1 

Au;«.-a;o^, 

xVrr 

where ojo ~ 3657 cm“^ is the vibration frequency for 
an isolated water molecule. The relative increase of 
the elastic constant AKrrV^rr^ at roi 02 = 2.8 A 
amounts to —0.05, i.e. Auj m —91.43 cm“^. Hence, 
the frequency shift sign for a dimer correlates with 
those in liquid water and ice. The shear moduli are 
identical by order of magnitude, but, nevertheless, 
they are considerably different. This fact has a sim¬ 
ple qualitative interpretation. The total electric field 
that acts on a water molecule in the liquid is, on 
the average, a little larger than that acting from the 


neighbor molecule in the dimer. An insignihcant in¬ 
crease of the electric field strength in the liquid is 
connected with a weakly ordered arrangement of the 
centers of mass of molecules and the orientations of 
its nearest neighbors. As a consequence, owing to the 
superposition principle, only a weak strengthening of 
the electric field in the molecule occurs. The opposite 
situation takes place in ice. 

Let us discuss the change of elastic constant in the 
standard dimers (see Fig. 5) at = 2.85 A. This 
is a distance between the oxygen atoms of water 
molecules in an argon matrix. According to our cal¬ 
culations, the relative increase of the elastic constant 
AK^rVK^r^ for this dimer configuration amounts to 
—0.0443 for the fixed orientation of molecules as in 
the standard dimer and to —0.0432 if the orienta¬ 
tions of molecules can be adjusted. These variations 
of the elastic constant correspond to a valence vi¬ 
bration frequency of 3576 cm“^ in the former case 
and 3578 cm“^ in the latter one. It should be noted 
that the relatively small value of orientational con¬ 
tribution is explained by the fact that the parame¬ 
ter of the system roi02 = 2.85 A is in the inter¬ 
val, where the dependence of the repulsion energy 
between molecules monotonically decreases. The fre¬ 
quency obtained from experiments in the argon ma¬ 
trix [38| equals 3574 cm“^. Therefore, we believe that 
it is possible to talk about a complete coincidence of 
calculated and experimental results. From our view¬ 
point, this is a powerful argument in favor that the 
frequency shift of valence vibrations has the electro¬ 
static origin. 

It is very important that the explanation of the 
frequency shifts different by magnitude is principally 
based on the superposition principle, the application 
of which to sharply directed and saturated irreducible 
hydrogen bonds is impossible. In order to make the 
substantiation of this fact more complete, we intend 


Table 6. Frequencies of symmetric 

valence vibrations of hydrogen atoms 

in the water molecule in vapor, water, and ice 



Uv , cm ^ 

Uv /, cm ^ 

I'ice, cm 1 

Experim. [151 

3657 

3490 

3200 

Experim. [32 

3656.7 

3280 


MST-FP [371 

3656 

3251 


SPC-FP 


3875 
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Fig. 7. Dependence of the ratio /fiQ on the distance 

between oxygen atoms in two water molecules 


to consider the frequency shifts of valence vibrations 
in ice and liquid water elsewhere. 

Not less demonstrative is the circumstance that one 
should expect a positive sigir of the frequency shift iir 
rarefied vapor, which directly follows from the behav¬ 
ior of AK^rVKrr^ in Fig. 6. This fact is also supported 
qualitatively by the experimental data on IR absorp¬ 
tion in rather rarefied water vapor . 

6. Influence of Neighbor Molecules 

on the Dipole Moment of a Water Molecule 

The dipole moment of an isolated water molecule is 
determiired as a sum of two oppositely directed vec¬ 
tors of dipole moments, p = /xh + Mo- The dipole 
moment /xh is determined by the spatial distribu¬ 
tion of the centers of oxygen’s negative charge and 
hydrogens’ positive charges, /xh = (;H(ri + r 2 ). The 
absolute value of dipole moment /xjj equals /xh = 
= 2qH?'OH cos (0/2) = 5.6281 D. The dipole momen¬ 
tum of the oxygen atom /xq emerges owing to the 
polarization of the oxygen anion’s electron shell in 
the electric field created by hydroren atoms in the 
water molecule. According to work [9|, it equals 

/xo = -aqn [1 -K{ri)] + [1 - K{r 2 )^ . 

It is easy to make sure that po = —3.7752 D. To¬ 
gether with /Xh, the following value is obtained for the 
absolute value of dipole moment /x: /x = /xh + /xo = 
= 1.8528 D. It completely agrees with the absolute 
value of dipole moment in an isolated water molecule. 

The variation of the dipole moment under the ac¬ 
tion of a neighbor molecule is one of the simplest man¬ 
ifestations of many-particle effects in the system. To 


estimate the influence of the second molecule, let us 
calculate the ratio /Xq Vmo between the dipole mo¬ 
ments of the oxygen atom in the pair approximation 
(/Xq^\ see work i) and in the isolated water molecule 
(/xo). The calculated dependence of /Xq Vmo on the 
distance between the oxygen atoms in two neighbor 
molecules is depicted in Fig. 7. 

One can see from Fig. 7 that the variation of the 
dipole moment of the oxygen atom under the in¬ 
fluence of the electric field created by the neighbor 
molecule does not exceed 1.5%. The same can also 
be said about the component /xh- From whence, it 
follows that the effects of electron shell overlapping, 
which are responsible for the change of the dipole mo¬ 
ment of a water molecule, are insignificant. This re¬ 
sult completely agrees with the conclusions of works 
Consequently, this means that the irreducible 
components of the interaction between molecules, 
which have to be regarded as hydrogen bonds, are 
much smaller in comparison with the energy of elec¬ 
trostatic interaction between molecules. 

7. Arguments in Favor 

of the Existence of a Hydrogen Bond 

In the previous sections, we presented the facts, 
whose explanation does not need the hypothesis 
about the existence of hydrogen bonds in water and 
other classical liquids according to L.A. Bulavin’s 
classification [l^. However, in this section, we de¬ 
scribe a phenomenon, which cannot be explained 
without attracting the hydrogen bond concept. This 
is the temperature dependence of the water heat ca¬ 
pacity. For convenience, this parameter will be reck¬ 
oned in the dimensionless units iq = 2Cv/ (A:bAa), 
where Cy is the heat capacity of a gram molecule, and 
Na Avogadro’s constant. These dimensionless units 
for the heat capacity will be called the number of 
thermal degrees of freedom. 

The latter differs from the standard number of de¬ 
grees of freedom in that the number of vibrational 
degrees of freedom in it doubles. Argon can serve as 
the simplest example here. In rarefied vapor, an ar¬ 
gon atom is described by three independent coordi¬ 
nates, which give its spatial position. The correspond¬ 
ing value of Xq also equals 3. In the crystalline state 
of argon, the number of standard degrees of freedom 
for its atoms also equals 3. However, iq = 6, because 
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each degree of freedom corresponds to the vibrational 
motion. 

The number iq of thermal degrees of freedom per 
water molecule consists of three components, 


.(tr) .(or) .(u) 

— *Q + *Q + ’ 


which correspond to the translational motion of the 
molecules, their rotation, and probable vibrations of 
irreducible hydrogen bonds, which are formed in all 
water phases 4dl4l| . 


From Fig. 8, one can see that the maximum value 
of heat capacity equals 6 for Ar in the liquid phase, 
12 for hydrogen sulfide, and reaches 20 and even more 
for water. The maximum value corresponds to a situ¬ 
ation where every of the ordinary degrees of freedom 
has a vibrational character. 

The heat capacity of liquid argon turns out some 
lower than in the solid state; nevertheless, it is close 
to 6. In the case of hydrogen sulfide, three orienta¬ 
tional degrees of freedom have also to be taken into 
account. If they have been vibrational, the maximum 
jQ-value would have been close to 12. Water would 
also have had this value for the number of thermal 
degrees of freedom per molecule if its molecules have 
not been bound with one another by means of hy¬ 
drogen bonds. However, actually, as one can see, the 
number jq for water exceeds this value approxima¬ 
tely by 6. 

From the physical viewpoint, this difference has a 
natural explanation: there are weak hydrogen bonds 
between molecules, which practically do not inter¬ 
fere with the rotational motion of the molecules, 
but make additional contributions to the water heat 
capacity. These contributions result from two trans¬ 
verse and one longitudinal vibrations of the hydrogen 
bond. Since every vibration corresponds to two ther¬ 
mal degrees of freedom, the excitation of all vibra¬ 
tions for only one hydrogen bond results in the heat 
capacity growth by 6. 


A more detailed analysis of the problem shows [17 1 
that, in order to attain a complete agreement with 
experimental data, it suffices to admit that every 
molecule of liquid water forms 2.5 hydrogen bonds 
near the triple point temperature and only one hy¬ 
drogen bond in a vicinity of the critical point. Those 
estimations are in quit e satisfactory agreement with 
the results of works |42l . l43l | , as well as with the results 
of computer calculations. 


8. Hydrogen Bond 

from the Viewpoint 

of the Chemical Bond Theory 

The hydrogen bond concept seems to appear for the 
first time in the work by A.R. Hantzsch in 1909 
[3|. Using the water molecule as an example, the 
sense of this new concept can be interpreted in the 
following way. The hydrogen bond is a new type of in¬ 
teraction between two water molecules. It acts along 
the O-H-0 line and is associated with the emer¬ 
gence of a specific interaction between the indicated 
groups of molecules that arises at certain distances 
between them. The introduced interaction is much 
stronger than the van der Waals one, but, at the same 
time, is much weaker than the covalent bond and the 
ionic interaction. In work [d^, an attempt was made 
to identify this specific interaction with the covalent 
bond between the hydrogen atom, which can hold two 
electron pairs about itself, and two electro neg ative 
atoms. L. Pauling criticized this approach |46l | and 
presented arguments in favor of the ionic nature of a 
hydrogen bond. Supposing that the hydrogen bond 
has a sharply directed character and is saturated, 
i.e. it does not agree with the superposition princi¬ 
ple, L. Pauling calculated the residual entropy of ice 
(= fee ln(3/2)) and showed that this value agrees well 
with experimental data (see also work 0). This re¬ 
sult promoted the further propagation of the hydro¬ 
gen bond concept while describing the properties of 
ice, water, alcohols, and so forth. 



0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 IJ 1.4 

t=T/X: 

Fig. 8. Temperature dependences of ieff for H 2 O, D 2 O, 
H 2 S, and argon systems in the liquid and vapor states at the 
curves of their coexistence. Experimental data are taken from 
work fjiil 


ISSN 2071-0186. Ukr. J. Phys. 2015. Vol. 60, No. 9 


969 


































I.V. Zhyganiuk, M.P. Malomuzh 


9. Discussion of the Results Obtained 

Let us summarize the results presented above, being 
based on the general ideas concerning the structure 
of intermolecular potentials in classical systems. We 
proceed from the fact that the simplest structure 
of the interaction between molecules is inherent to 
atomic systems of the argon type. The correspond¬ 
ing potential of the intermolecular interaction $(r) is 
a sum of an attractive component which is 

associated with dispersion forces, and a component 
<l?rep(?') describing the repulsion: 

$(r) = $dis(?’) + $rep(?’)- (14) 

Note that the known Lennard-Jones potential has just 
this structure. 

The systems consisting of molecules like N 2 are 
characterized by the loss of spherical symmetry, 
which is accompanied by the appearance of the appre¬ 
ciable an gula r dependence in the particle-to-particle 
potential |48|: 

$(r) $(r, n) = $dis(r', fl) -f 4’rep(?’, ^)- (15) 

This structure must also be inherent to the interac¬ 
tion potentials between water and alcohol molecules 



Fig. 9. Partial contributions to the interaction potential 
between water molecules in the linear dimer 



Fig. 10. Hydrogen bond potential in the linear dimer of 
water molecules 
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if the distance between them considerably exceeds the 
sum of their molecular radii. 

However, as the molecules get closer, their electron 
shells overlap, and a qualitatively new component 
of the intermolecular interaction, <i>H(r, H), appears, 
which is conventionally called the hydrogen bond en¬ 
ergy. In this case, 

$(r, ft) = $dis(L + ^repir, H) -f 

-f $£;(r,H)-f $H(Lf^)- (16) 

At the same time, the hydrogen bond energy, as a 
rule (see work (i^), is associated with the sum of two 
last terms, 

AH(r,H) = $£;(r,L!) + $H(r-,H), (17) 

and M.D. Sokolov was the first who attracted atten¬ 
tion to this fact for the first time [ij. 

In order to estimate the irreducible contribution 
to the hydrogen bond ener gy, the following ap¬ 
proach was proposed in work |49l |. From the conven¬ 
tional viewpoint, the ground state energy of a wa¬ 
ter dimer is determined by the “hydrogen bond” en¬ 
ergy Ep[(rd, Hd), where the subscript d indicates that 
the distance between the oxygen atoms of two wa¬ 
ter molecules and the angle values correspond to the 
dimer configuration. On the other hand, the prop¬ 
erties of dimers are described quite well, by us¬ 
ing phenomenological intermolecular potentials of the 
SPG 0, SPC/E [EH, TIPS [nl, SD 0, GSD 0 , and 
so forth types. 

The dimer configuration is determined from the 
condition of interaction energy minimum. Proceeding 
from this fact and expression m, the energy of an 
irreducible hydrogen bond can be estimated with the 
help of the relation 

$H(r, H) = $(r, H) - 

- [4>dis(r', H) -I- $rep(L H) -I- $M(r’, H)], (18) 

where $m(?', U) is the multipole approximation to the 
electrostatic interaction energy. The multipole mo¬ 
ments of a water molecule are assumed to be de¬ 
termined independently of effective charges govern¬ 
ing the behavior of the phenomenological potential 
<i>(r, fl). The magnitudes of effective charges, as well 
as some other parameters of phenomenological poten¬ 
tials, are determined by fitting the model to dimer 
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parameters, which are determined experimentally or 
within quantum-chemical methods. Hence, it is the 
difference on the right-hand side of Eq. dT51) that de¬ 
scribes the component emerging owing to the over¬ 
lapping of electron shells in water molecules. In work 
|49| . the intermolecular potential <i)(r, H) was sim¬ 
ulated by the generalized Stillinger-David potential 
^’GSD(r'), so that 

$H(r, H) = $gsd(?', H) - 

- [$rep(?’, H) -f H) -I- H)]. (19) 

The general behavior of the potentials 4 >gsd(?’5H), 
^M{r, H), and $<iis('r, H) for the molecular orientation 
typical of the dimer are exhibited in Fig. 9. 

The hydrogen bond potential for the same config¬ 
uration of water molecules is shown in Fig. 10. This 
is a short-range potential that emerges owing to the 
overlapping of electron shells and has a quantum- 
mechanical origin. It is this potential that should be 
interpreted as the hydrogen bond potential in wa¬ 
ter. Its depth, by order of magnitude, is identical to 
the depth of the dispersion force interaction poten¬ 
tial between water molecules, but it is substantially 
smaller by magnitude than the multipole interaction 
potential. As a result, the contribution made by hy¬ 
drogen bonds to the thermodynamic potentials of wa¬ 
ter can be taken into account in the framework of the 
thermodynamic perturbation theory. At the qualita¬ 
tive level, this circumstance completely agrees with 
the similarity of the thermodynamic functions of wa¬ 
ter and argon at their coexistence curves. 

The presented estimations of the magnitudes of the 
electrostatic interaction and hydrogen bond contribu¬ 
tions to the energy of intermolecular interaction were 
confirmed many times in works . In particular, 
it was shown in work that, by order of magnitude, 

g2)| ^ W 
|t^H(gl,g2)| 2 

This means that the influence of exactly hydrogen 
bonds is taken into account in a natural way by the 
thermodynamic perturbation theory (T^ . [2^ . 

While analyzing the thermodynamic properties of 
liquids and solutions with the use of statistical theory, 
it is necessary to take into consideration the fact that 
molecules permanently rotate and, at the same time, 
hinder their rotational motion. The characteristic pe¬ 
riod of the thermal rotational motion of molecules 


turns out much shorter than the characteristic time 
of changing the configurations formed by the transla¬ 
tional degrees of freedom. Therefore, the thermody¬ 
namic properties of liquids are mainly governed by the 
potentials UAiri 2 ) averaged over all angular variables 
(here, ri 2 is the distance between the centers of mass 
of molecules). In works [i3;[23j it was demonstrated 
that such an averaged potential has a structure of the 
Sutherland potential. 


UA{ri2) => 


foo, ri 2 < ro 

\Ua{ri2), ri2 > ro. 


( 20 ) 


where Uair) is the attraction potential, which de¬ 
creases at large enough distances following the law 
1/^12- With the same accuracy, the averaged poten¬ 
tial can be approximated by the Lennard-Jones po¬ 
tential 


UA{ri2) => C/Lj(ri2), 


C^Lj(ri2) = -4e 


a2 


-) -V 

ri2y Vri2 


6n 


( 21 ) 


The thermodynamic properties of liquids consist¬ 
ing of anisotropic molecules, owing to the rotational 
motion of those molecules, are similar to the corre- 
spon ding p roperties of atomic liquids of the argon 
type 17|,|23| . Small differences between those proper¬ 
ties stem from weak angular correlations, which can 
be taken into account with the help of perturbation 
theory. The most noticeable angular correlations are 
observed in the temperature interval of the super¬ 
cooled liquid state and at the formation of an instant 
local structure in liquids. This circumstance is a pre¬ 
condition for the emergence of special points in aque¬ 
ous alcohol solutions. 

From the qualitative viewpoint, such a change of 
priorities is not justified, because the analytic con¬ 
tinuation of the components $(iis('r, H) and <I>m(?') H) 
into the region, where electron shells overlap is not 
accompanied by the appearance of effects that would 
violate the requirement of continuity in the potential 
behavior. For this reason, it is desirable that the hy¬ 
drogen bond potential should be defined in another 
way. According to the continuity requirement, the hy¬ 
drogen bond potential will be defined, by using the 
formula 


4>(r, H) = $dis(r, H)-|- 

-|-$r.(r, H) -I- $m(?', H) -I- $H(r, H). 
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In the region where the electron shells overlap, the 
functions <i)dis(?', and fl) should be sub¬ 

stituted by their corresponding continuations from 
the region, where their application does not invoke 
doubts. 

The formation of dimers and multimers of higher 
orders in vaporous and liquid water is one of the 
most characteristic manifestations of hydrogen bonds, 
which are responsible for the specificity of the in¬ 
teraction between molecules. In other words, the 
study of the properties of dimers in water provides 
us with direct information concerning the proper¬ 
ties of hydrogen bonds. This circumstance gives us 
an exact instruction on how one must approach 
the research of the properties of the interaction be¬ 
tween molecules in water and the hydrogen bond 
formation itself. The main stages of such an ap¬ 
proach are as follows. First, to describe the en¬ 
ergy of interaction between two water molecules, 
the most suitable phenomenological model potential 
is selected among those that describe the ground- 
state energy of a dimer the most successfully. At 
the second stage, the obtained dependence of the 
interaction energy of water molecules on the dis¬ 
tance between them is compared with the interac¬ 
tion energy determined from the asymptotic mul¬ 
tipole series expansion. Finally, in order to deter¬ 
mine the dependence of the hydrogen bond energy 
on the distance between molecules, the difference be¬ 
tween the energies of the model potential and the 
sum of the dispersion and multipole components 
is calculated. This difference is expected to be dif¬ 
ferent from zero only within a certain vicinity of 
the equilibrium distance between water molecules 
in the dimer. 

In order to describe the interaction between 
molecules in a dimer, we use the generalized Stil- 
linger-David potential proposed in work Q. It is a 
soft potential, whose parameters can be changed due 
to the interaction with neighbor molecules. This is a 
very important circumstance, which cannot be taken 
into account in the majority of phenomenological 
model potentials [IlJ3- Unlike the original Stillin- 
ger-David potential 9|, its generalized variant (GSD) 
more adequately involves the behavior of screening 
functions, which describe the effects of electron shell 
overlapping. In addition, the Stillinger-David poten¬ 
tial was corrected with respect to it asymptotic be¬ 
havior at large enough distances between molecules, 


when it should be determined by the dipole-dipole 
interaction. 

One can see that the depth of the irreducible com¬ 
ponent d>H(r, U) in the interaction potential between 
water molecules resulting from the overlapping of 
their electron shells does not exceed {2^3)k^T^. By 
order of magnitude, it is close to the dispersion 
component $(jis(r, 0), but is substantially smaller in 
comparison with the multipole interaction potential 
$M(r,U) Ri (7-8)fcBT™. 

Actually, the insignificant depth of the potential 
well formed by the hydrogen bond results in that the 
contributions of the interaction potential component 
to the thermodynamic potentials and the kinetic co¬ 
efficients can be taken into account in the framework 
of perturbation theory. In addition, the temperature 
behavior of main thermodynamic parameters of wa¬ 
ter such as the fractional molecular volume, evap¬ 
oration heat, and others, has an argon-like charac¬ 
ter with a quite satisfactory accuracy. These conclu¬ 
sions are completely confirmed by the results of works 


The hydrogen bond potential plotted in Fig. 10 cor¬ 
responds to the relative orientation of water molecules 
in the equilibrium dimer configuration only. In prin¬ 
ciple, there are no complications for the construction 
of the potential U) at all other relative orien¬ 

tations of water molecules, since the angular depen¬ 
dences for the potentials $gsd(?', U) and U) are 

known for arbitrary angles. 

The conclusion about a weak deformation of elec¬ 
tron shells and, as a consequence, the formation of 
weak irreducible hydrogen bonds is obviously sup¬ 
ported by the results of work [Hi , in which a redistri¬ 
bution of the electron density was analyzed, by using 
the methods of scanning tunnel microscopy. 

It should be noted that the thermodynamic prop¬ 
erties of water are determined by the potentials av¬ 
eraged over the angles, which is a consequence of 
the rotational motions of water molecules. Owing 
to this averaging, as was shown in work [H|, the 
potential well depth of a hydrogen bond addition¬ 
ally decreases, which gives rise to a correction of 
the argon-like dependences for the thermodynamic 
parameters, whose relative magnitude does not ex¬ 
ceed 5% [l3, [H|- At the same time, the hydro¬ 
gen bonds manifest themselves directly in the wa¬ 
ter heat capacity [13 . Another important circum¬ 
stance falling beyond the scope of our consideration 
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is the adequate account for the environment influ¬ 
ence on the character of the hydrogen bond poten¬ 
tial. We are planning to consider this issue in detail 
elsewhere. 

This work would be impossible without regular con¬ 
sultations over the years with Profs. T. V. Lokotosh, 
G.G. Malenkov, Yu.I. Naberukhin, G.O. Puchkovska, 
and V.E. Pogorelov. We are also grateful to our coau¬ 
thors P. V. Makhlaichuk and S. V. Lishchuk. While 
carrying out this and other works in this direction, we 
felt the constant support of Academician of the NAS 
of Ukraine L.A. Bulavin, for which we are sincerely 
thankful to him. 

The results of those works were reported at various 
seminars and conferences. After the first reports were 
made, a careful criticism of attendees was gradually 
changed in favor of our viewpoint. 
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Translated from Ukrainian by O.I. Voitenko 


LB. IKusanwK, M.II. MaAOMyofc 

OI3HHHA nPHPO^A BO^HEBOPO SB’HSKV 

P e 3 lo M e 

y po6oTi flOCjii,i];acyeTi>CH (J)i3HHHa npnpofla xa KopeKTHicTb 
oaHaueHHH Bo;i;HeBHx 3B’B:3KiB. AnajiiayeTBCH, nepm 3a Bce, 
BnjiHB ocTannix na noBe;i;iHKy TepMo;i;HHaMiHHHx, KineTH- 
HHHX Ta cneKxpocKoniuHHx BjiacTHBOCxen bo^h. HoKaaaHO, n],o 
cnpHHHHTTH BO^^HeBHx 3B’B:3KiB HK TOCTpoHanpaBjieHHx xa na- 
CHueHHx MicxKiB, HKi BHHHKaioxB Niixc MOjieKyjiaMH bo;i;h, e 
necyMicHHM 3 noBe;i;iHKOio cneLi,H4)iHHoro o6’eMy xa xenjioxn 
BHnapoByBaHHH, a xaKo:>K Koec^ipienxiB caMO^H4)y3ii xa Kine- 
MaxHHHoi" 3cyBHOi B’H3KOCxi. Ha ;i;o;i;axoK ;],o pnoro noKaaano, 
u],o 3MiHa pnnojiBHoro MOMenxy MOjiexyji Bopn, a xaKo:ac 3cyB 
uacxoxH BajienxHHx KOjinBanb rippoKCHjiBHoi" rpynn noBnicxio 
noHCHioioxbCB: na ocHOBi ynBjienb npo ejieKxpocxaxHuny npn- 
popy BopneBoro 3B’B:3Ky. Pa30M 3 xhm, xeMnepaxypni 3ajie:acHO- 
cxi xenjioeMHOCxi Bopn xa li* napn uixKO BKaayioxb na icnyBan- 
nn cjiabKHx Bopneanx 3B’H3KiB. AnajiiayiOHH BjiacxHBOCxi pn- 
Mepy Bopn, noKaaano, mo bhccok cjiabKHx Bopneanx 3B’H3KiB 
y enepriio ochobhofo cxany pniuepy e npH6jiH3HO a 4—5 paaia 
MeHmHM y nopiBHHHHi 3 eneprieio ejieKxpocxaxHUHoi" B3aeMop;ii 
MiTK MOjieKyjiaMH Bopn. HipcyMOByioun pe3yjibxaxH, pobuxacn 
BHCHOBOK, mo xBKy caiuy npnpopy Maioxb BopHeai 3b’h3kh b 
ycix inmnx annapKax, pe bohh BHHHKaioxb. 
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